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We consider a magnetic superconductor (MS) with a spiral magnetic structure. On the basis 
of generahzed Eilenberger and Usadel equations we show that near the boundary of the MS with 
an insulator or vacuum the condensate (Gor'kov's) Green's functions are disturbed by boundary 
conditions and differ essentially from their values in the bulk. Corrections to the bulk quasiclassical 
Green's functions oscillate with the period of the magnetic spiral, 2n/Q, and decay inside the 
superconductor over a length of the order v/2-kT (ballistic limit) or yj D/tvT (diffusive limit). We 
calculate the dc Josephson current in an MS/I/MS tunnel junction and show that the critical 
Josephson current differs substantially from that obtained with the help of the tunnel Hamiltonian 
method and bulk Green's functions. 

PACS numbers: 74.50.+r.De, 74.45. +c 



I. INTRODUCTION 

It is known that in some compounds the superconducting order can coexist with a magnetic order of the ferromag- 
netic or antiferromagnetic type. For example, in ternary rare-earth compounds such as (RE)Rh4B4 and (RE)Mo6X8 
(X=S,Se) the superconducting and magnetic ordering coexists in a narrow temperature range (see the review ^1] and 
a more recent paper Q and references therein). In ErRh4B4 superconductivity takes place in the interval 0.7 < T 
< 0.8 K, and the magnetic ordering arises below rm=0.8 K. In HoMoeSg the magntic ordering occurs below Tm— 
0.74 K, whereas superconductivty exists in the temperature range 0.7 < T < 1.8 K. Besides, the superconducting and 
magnetic order is realized in the layered perovskite ruthenocuprate compound RuSr2GdCu208 [3,13, H, 0, 0] • this 
compound an antiferromagnetic order and, perhaps, a weak ferromagnetism take place. 

A uniform magnetization is impossible in a bulk superconductor as the magnetic field destroys superconductivity. 
In order to explain the coexistance of ferromagnetism and superconductivity, Ginsburg and later Anderson and Suhl 
supposed that this coexistance is possible in case of a domain or spiral magnetic structure [1, 0]. The period of 
the magnetic structure has been calculated in Ref.0] (see also Ref.[l2|), and on the order of magnitude it is equal 
to Im ~ 2n{^Qkpy/^ /kp, where kp is the Fermi momentum and = ^^f/ttAq is the correlation length in a clean 
superconductor. For example, in HoMoeSs the wave vector of the periodic magnetic structure Q « 0.03 [^[Tol.[T]|. 

As is well known, many characteristics of a superconductor (the critical temperature, the density-of-states etc) 
can be calculated if the Green's functions of the system, including the anomalous ones (or Gor'kov's functions), F, 
are found [l^. These functions for a magnetic superconductor (MS) with a spiral structure have been obtained in 
Ref. [l3| • In this case the functions F depend on the center-of-mass coordinate and momentum direction so that the 
system is anisotropic. Long ago, it was established that surface effects are essential for finite anisotropic sam ples 
such as anisotropic superconductors and high Tc superconductors with d wave pairing (see, for example, Ref.|14l| 
and also the review |15| and references therein). In particular, the order parameter may be suppressed near the 
superconductor/ vacuum or superconductor/insulator {S/V or S/I) interface. A high impurity concentration leads to 
averaging the Green's functions in the momentum space so that in the diffusive limit, characteristics of the system 
do not depend on the sample size. 

In this paper, we show that the surface effects are important in MSs with a spiral magnetic structure. In particular, 
the Green's functions of the system are disturbed by boundary conditions at the S/V or S/I interface in samples with 
any impurity concentration. Corrections to the bulk Greens functions due to boundary conditions oscillate in space 
with the period 2tt/Q and decay from the interface over a length of the order ~ v/2ttT in the ballistic limit and 
of the order — \JD/ttT in the diffusive Hmit. 

The surface effects become very important in the cases when one needs to know the Green's functions near the 
interfaces. For instance, the Josephson current Ij in an MS/I/MS tunnel junction is determined by the values of the 
Green's functions near the MS/I interface (/ stands for an insulating layer). The Josephson current Ij in the MS/I/MS 
junction with a spiral magnetic structure was calculated in Ref. jl6| on the basis of the tunnel Hamiltonian method. 
The authors used the Gor'kov's functions calculated in Ref.[l2| for an infinite MS with a spiral magnetic structure 
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in the ballistic limit. They have obtained that the Josephson critical current Ic depends on the angle 6 between the 
magnetization directions in both MSs near the interface and calculated the dependence of 1^ on different parameters 
of the junction (the exchange field, the wave vector of the spiral, Q, etc). It has been established that at some values 
of parameters the critical current becomes negative (tt - state). We will show here that, although the current Ic 
indeed depends on in a way similar to that in Ref. J^j , the dependence of Ic on various parameters is completely 
different. The point is that the tunnel Hamiltonian method is not applicable to inhomogcneous superconductors 
and, in particular, to M^s with a spiral magnetization. In order to calculate /c, one has to solve the Eilenberger 
or Usadel equation with boundary conditions at the MS/I interface. It turns out that the Green's functions at the 
MS/I interface differ essentially from their values in the bulk, and correspondingly the Josephson current also differs 
substantially from its value obtained on the basis of the bulk Green's functions. 

The structure of the paper is the following. In Sec. II, we analyze the ballistic case. Using the Eilenberger equation 
generalized for the case of the MS with a magnetic spiral, we find the spatial dependence of corrections to the bulk 
Green's functions. In Sec. Ill, the diffusive case will be considered. Using a generahzed Usadel equation complemented 
by boundary conditions at the MS/I interface, we calculate the Josephson current in MS/I/MS tunnel junction and 
compare the obtained critical Josephson current Ic with that obtained on the basis of the tunnel Hamiltonian method. 
In Sec. IV, we discuss the obtained results. 



II. BALLISTIC CASE 



We consider a MS with a spiral magnetic structure. The exchange field acting on free electrons is assumed to 
lay in the (y, z) plane and to rotate in space with the wave vector Q; that is, the vector of the exchange field 
is: h — ft,(0, sina(a;), cos Q!(x)) with a = Qx + 9, x > 0, {0 is the angle between the magnetization and z-axis at 
X — 0). The superconducting order parameter A is taken into account in the mean field approximation: A — 
As X]p(^T,pV'i,-p)j the singlet pairing is assumed. The Eilenberger equation is derived in a standard way (see, 
for example, [H [H, H E!, HI ) . The main difference between the cases of an ordinary, nonmagnetic superconductor 
and MS with a spiral structure is that the quasiclassical Green's function g in the latter case is a 4 x 4 matrix in the 
Gor'kov-Nambu and spin space. This equation has the form 



ivVg + cj[f3 (g)(5-o,g] +i[h{x)S,g] 



A( 



(z/2r) m,9]=0 



(1) 



where v is the Fermi velocity, S = (cti,5'2,T3 (giCTa), ak,fk are the Pauli matrices in the spin and Gor'kov-Nambu 
space, and (To,'?o are the unit matrices. The square and angle brackets mean the commutator and averaging over 
angles, respectively, and r is an elastic scattering time. In order to exclude the coordinate dependence of the third 
term in Eq.([T]), we perform a transformation (see Ref. [2(1] ) 



g = U®g„m+, (2) 

where U — tqiE) d'o cos(a/2) + i sin(Q!/2)f3 (g) ai is an operator corresponding to a rotation in the spin and particle-hole 
space, and gn is a new matrix. Then Eq.([T|) acquires the form 



vfid^g + [t3 (E> {ujo-q + iha3,g] + iv^i{Q/2)[T3 (g>ai,g] 



' <^3,9 



(z/2t) [{g),g]^0 



(3) 



where fi — Px/p- The subindex "n" is omitted. From the physical point of view, the transformation given by Eq.Q 
means the transition to a rotating coordinate system, in which the magnetization vector is directed along the z-axis. 
That is why the exchange field h in Eq.Q contains only the z-component. 

For simplicity, we restrict the consideration with the case of temperatures close to the critical one of the supercon- 
ducting transition, Tc. In this case the matrix Green's function g may be represented in the form 



g = signu) • T3 (g) ctq 4- / . (4) 

where the anomalous (Gor'kov's) matrix function, / , is assumed to be small, that is, all elements of this matrix are 
small. The first term is the normal, matrix Green's function in the Matsubara representation. 
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In this Section, we consider the ballistic case, i.e. we suppose that r ^ oo. Substituting the matrix g from Eq.Q 
into Eq.(IT]), we come to the equation for the anomalous function / 

vfiT3(E)dxf + ivfi{Q/2)[ai, f]+ + 2w/ + ih[a3, /]+ = f2 (8) asAsignu. (5) 
We represent the matrix / in the form 

/ =/®f2+F®fi (6) 
where / and F are matrices in the spin space that can be represented as a sum of Pauli matrices 

/^E/fc'^fc; ^^E^fc'^fc (7) 
fc fc 

where fc = 0, 1, 3. 

Eq. ([ni) is a system of linear equations with respect to coefficients fk and Fk . The solution of these equations consists 
of a part, fk and Fk, constant in space and a nonhomogeneous part, dfk{x) and 5Fk{x). The latter part arises if there 
are nontrivial boundary conditions in the problem. The homogeneous part is a solution for an infinite sample when 
boundary conditions can be ignored. The homogeneous solution can be easily found. It has the form 

j_ A(e|+cj2) ihAsignuj 



where eg = All other coefficients (i.e. fi,Fk) equal to zero. The coefficient /a is the amplitude of the 

singlet component, and the coefficient fo is the amplitude of the triplet component with zero projection of the total 
spin of a Cooper pair on the 2;-axis (in the rotating coordinate system), Sz = 0- The singlet component is an even 
function of uj, while the triplet component, /o, is an odd function of lo 20] . One can see that the exchange field, h, 
suppresses the amplitude /a, whereas at a sufficiently large wave vector of the spiral Q, the amplitude /a is restored 
to the value A/|ct;| which is the amplitude of the condensate function in a nonmagnetic superconductor. Note that 
the authors of Ref.[16,] used only bulk solutions in the laboratory coordinate frame. These functions may be reduced 
to the quasiclassical Green's functions in Eq.®. 

The function /a determines a change of the critical temperature of the superconducting transition, T^, due to the 
exchange field h and wave vector of the magnetic spiral Q (see, for example, the review articles [l9l[20|) 



TrO — Tr 



(9) 



where Tco is the critical temperature in the absence of the exchange field h. It is seen that with decreasing the spiral 
period, 2tt/Q, the suppression of the critical temperature is reduced and at vQ » h the critical temperature is the 
same as in a nonmagnetic superconductor, i.e. Tc —^ Tco- 

Now we turn to the calculation of corrections dfk and SFk that arise due to boundary conditions and depend on x. 
Note that if the correction Sfs is not small compared to /a, a correction to the order parameter SA(x) will not be small 
as well. This circumstance makes the problem rather complicated because Eq. ([5]) becomes a system of six equations 
with the right-hand side which depends on x. In order to simplify the problem, we assume that the correction (5 /a is 
small and we can neglect a variation of A in space. We will see below that in a general case (5 /a may be comparable 
with /a. In this case our results are correct up to a numerical factor of the order unity. In the next Section, we discuss 
the validity of the obtained results in more detail. 

Thus, in order to find the corrections Sfk and 6Fk, we have to solve a system of homogeneous linear equations ([5]) 
without the right-hand side. Substituting the expansions ([7]) with Sfk and SFk as the coefficients of these expansions 
into Eq.® with 6 A = and representing the coordinate dependence of these coefficients in the form {Sfk,SFk} ^ 
exp{Kx), we obtain a system of six linear equations. One can see from these equations that the coefficients /i and 
Fo.a are antisymmetric functions of whereas the coefficients /o,a and Fi are symmetric functions of /j,. We do not 
write down these equations as they are rather cumbersome. Instead of this, we write the determinant of the system 
which determines the eigenvalues Kj. It is reduced to a cubic algebraic equation 
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A^(l - z)[A^(l - A^) + z^] + (1 - A^)(l + zf{\^ - z) = (10) 

where z = e^/f^^, Ah = h'^/n'^, Xq = e^/fi^, Q'^ ^ uj^ + h'^ + e|, = fivK/2,and eg = fivQ/2. 

In order to find the eigenvalues, one has to solve this equation. We consider the most interesting case of large 
energy eg : eg >> T, h. In this case the critical temperature Tc is close to Tco- The solutions of Eq. fTO)) are 



Zl 



= A,^, Ki = 2\uj\/{v\n\) 



(11) 



and 



Z2,3 = -1 ± K2,3 = ±«(3 - 2\Lu\/{v\fl\) (12) 

Therefore, the eigenfunctions corresponding to K2,3 oscillate in space with the period of the spiral and decay over 
the distance of the order S.T = V/2TTT. The eigenf unction, which corresponds to ni, decreases monotonously from the 
interface over the correlation length ^t- 

The amplitudes fk and Fj. may be found from boundary conditions at the MS/V or MS/I interfaces [s^l 



/ (m) - / (-M) = (13) 

which read that the antisymmetric part of the Green's function should turn to zero at the MS/I interface. One can 
solve the corresponding equations and find the amplitudes fk and Fk. However we will not do that for two reasons. 
First, the corresponding expressions are cumbersome. The second and more important reason is that the surface 
effects are displayed near the interface at which a random (diffusive) scattering takes place. Therefore, the ballistic 
case considered in this Section is not relevant to this situation. In the next Section we consider a more realistic case 
of a sample with a high impurity concentration (dirty case). We will find the eigenvalues Ki and the amplitudes of 
eigenfunctions. One can show that the structure and form of the dependencies of the functions fk and on fi and lo 
are qualitatively the same in both cases, ballistic and diffusive. The only difference is that, whereas in the diffusive 
case only the zero and first terms in the expansion in spherical harmonics are important, in the ballistic case the 
dependence on is more complicated. 



III. DIFFUSIVE CASE 



In this Section, we consider the influence of the boundary on the condensate functions assuming that the impurity 
concentration is high and the condition / << 2'k/Q,£^t is satisfied, where / = ut is the mean free path. In this case 
the part of the condensate function fasm antisymmetric in the momentum space is expressed through the symmetric 
part via the well known expression [l^, [l^, [2^ l2ll . [23 | 



f asm = {Px/\p\)T3sgnuj{d^f +i(Q/2)f3[(Ti,/]+) (14) 

where r^sgnuj is the ordinary quasiclassical Green's function in the normal state (see Eq.Q). The second term arises 
as a result of the transformation the term [(Ti,/]+ means anti-commutator. One can see that the asymmetric 
part has the opposite parity in to compared to the symmetric part f ; if f is an odd function of uj, then / asm is an 
even function of uj and vice versa [20| . In the simplest case of ordinary BCS superconductors the symmetric function 
near Tc is equal to / = fa iJoA/|u;|, i.e. is an even function of uj. Obviously the antisymmetric part / asm is an odd 
function of uj. This issue is discussed in detail in Refs.psj. 

We assume again that the temperature is close to Tc. The symmetric part of the condensate function / after the 
transformation Eq.([2]) obeys the equation [20l | 



D{-dlJ+^if + ai®f®ai) + i^r^iai, dj] + } + 2\uj\f + z/i^[a3, /]+ = 2f2 ® agA. (15) 
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Here = sgnojh. As follows from Eqs. (|13ll4p the boundary condition has the form 

5,/ +z(Q/2)f3[ai,/]+ = (16) 

This means that the spiral axis is assumed to be perpendicular to the MS/V or MS/I interface. 

One can see that a coordinate-independent solution for Eq. (fT5|) satisfies the boundary condition only if Q = 0. If 
Q is not zero, the anti-commutator [(T3, /o(7o]+ 7^ 0, and therefore dxf also differs from zero at the boundary. 

We have to solve Eq. (fT5|) with the boundary condition p6| . The uniform solution again has the form ([8]) with 
eg = DQ^ /2. The correction Sf = f — T2(/3(T3 + /o^o) satisfies the uniform equation (jisp and may be represented 
in the form ^Eil, where only the coefficients /o,3 and Fi differ from zero, that is 

/ = (/3<3-3 + /o<To) «) fa + FiCTi (g) fi (17) 

We look for a solution in the form of exponentially decaying functions: 5f ^ exjp{Kx) with Hen < 0. The 
determinant of the system of Eqs. p3)l has the form 

[(1 + zf + 2A^(1 -z) + Xl]{X^ ~z) + Xl{l + - z) = (18) 

where z = [n/Qf, = 2\uj\/DQ'', Xh = 2h^/DQ^. 

Again we consider the most interesting case of small Xi^,h which seems to be relevant to the experiment pi)| : 
{^Lu, A/i} << 1, i.e. {T, h} « DQ^ ^2^. In this hmit the eigenvalues are 



zi=X^+ Xl; Ki = -Q^X^ + Xl (19) 

and 



Z2,3 = -1 ± i^2{2X^ + Xl)- ^2,3 = - Q^{2X^ + A^J/2 (20) 
Thus, the correction 5/3 may be written as 

^fi{x) = 0.1 exp(Kix) -I- a+ exp(K_|_x) -|- a_ exp(K_a;) (21) 



where k+ = K2 = +iQ - gV(2A„ + A^)/2 and K,^ — = K*^. The first term decreases monotonously inside the 
superconductor, whereas the second and third terms oscillate with the period 27: /Q and decay over the length of the 
order of min{^T, \/ iDQ^/h){D/h)}. The corrections, Sfo and Fi, have the form 

^foix) — —iXhtti exp{Kix) — {iXh)^^[{l + ia)a+ exp{K+x) -f (1 — ia)a- exp(K_x)] (22) 
Fi{x) = 2iXh^/ziai exp(Kia::) — A^"'"[(l + ia)a+ exp(K+a;) — (1 — ia)a^ exp(K_x)] (23) 



where a = ■y/2(2Aij 4- Xf/j. The coefficients oi and a± are found from the boundary condition ([T 

A^ + cty/zi/ Z 2 a 

Making use of Eqs. (|2m23)) . one can obtain the values of the condensate function at the interface /(O) that determine, 
for example, the Josephson current in MS/I/MS junction. We find 
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FIG. 1: Contributions of the singlet (solid line) and triplet Sz ~ (dotted line) components to the Josephson critical current 
as a function of the exchange field. This dependence has been obtained on the basis of bulk quasiclassical Green's functions 
and corresponds to the tunnel Hamiltonian method. The critical current Ic{h) and the exchange field h are plotted in units 
/c(0) and DQ^ /2, respectively. The temperature is chosen equal to T = O.IDQ^ /2'jt. 

and 

i^i(0) = ~z^/3(0) (26) 
a 

where the amplitude of the bulk singlet component can be expressed in terms of the parameters X^jji 

f + 1) 2A 

\^{\^ + l) + \lDQ^' '^"> 

In the considered limit, X^^h <C 1, the function /a is close to the value of the singlet component in an ordinary 
(nonmagnetic) superconductor. The exchange field, which tries to destroy Cooper pairs, is effectively averaged due 
to rotation of the magnetization vector. 

Now we discuss the co nditions under which the obtained results are valid. Consider first the case of a thick sample 
{d ^ ^GL — 1.2\AD/r(T/A), where d is the thickness of the sample and £,gl is the Ginsburg-Landau correlation 
length) One can see that if Xl » X^{n = 0), i.e. » {'kT)DQ'^, the value of /3(0) is /3(0) = hV^I{2 + \/2) w 
0.41/3, i.e. the singlet condensate function at the interface differs from the bulk value by a numerical factor of the 
order 1. In the limit X\ « X^j the singlet component is almost constant in space so that /3(0) « fs. Therefore 
our results are valid in this limit. However, our results are also correct if the thickness of the sample d is less than 
the Ginsburg-Landau correlation length. In this case, the order parameter A is constant in space [2^ so that our 
assumption about the coordinate-independent A is fulfilled and the obtained results are exact. 

Let us discuss the meaning of the component of the condensate function / . As we said above, the function /3(0) 
is the amplitude of the singlet component at the interface. The function /o(0) is the amplitude of the triplet 5^=0 
component. One can see that both functions, /o and Sfo, (the bulk value and the correction due to the surface effects) 



7 




FIG. 2: Contributions of the singlet (solid line) and triplet (dotted lines) components to the Josephson critical current as a 
function of the exchange field h. The upper and lower dotted curves correspond to the triplet \Sz\ = 1 and Sz = Q condensate 
components. These curves are plotted on the basis of Ea. (l28|l . The normalization units are the same as in Fig.l. 



are small compared to the singlet value in the considered limits, X^i.h « 1- The function Fi is the amplitude of the 
triplet component with \Sz \ = 1 in the rotating coordinate system. In the bulk, it is equal to zero. Just this component 
penetrates the ferromagnet over a long distance in S/F structures with a rotating magnetization [20, 26, 27, 28, 29]. 
This triplet component -Fi(O) is of the order of the singlet component in the bulk, /s, at << X\ and less than /a 
at > ^ . 

Knowing the quasicalssical Green's functions at the MS/I interface, we can calculate the dc Josephson current Ij 
in a MS/I/MS tunnel junction consisiting of two MSs. The Josephson current in this junction is expressed in terms 
of the components /o,3 and Fi at the interfaces MS/I, i.e. at a: = (see the Appendix) 



Ij = /,sin(p, h = (ei?B)-i(2^T) ^^(/f (0) + cos0[/2(O) + ^^(0)]) (28) 

a)=0 

where Rb is the resistance of the junction in the normal state, is the phase difference, uj — 7rT(2n + 1) is the Mat- 
subara frequency, and 6 is the angle between the magnetization vectors in the right and left magnetic superconductors 
at the interfaces. Since we are interested in the Josephson current in the lowest order in the parameter ii^^, the 
functions fo,3ix) and Fi(x) should obey the boundary conditions (jl6p that correspond to the limit Rb — *■ oo. These 
functions are given by Eqs. (|25ll25)) . 

A formula, which resembles Ea. ((28)) . was obtained in Ref. [16] on the basis of the tunnel Hamiltonian method. 
What is the difference between these two formulae? First, the term -Ff (0) is absent in Ref.[l^. Second, instead of 
terms /o3(0), in Rehfl^ there are terms /03 corresponding to the bulk solutions. This difference leads to essential 
consequences. In particular, the conclusion made in Ref. 1 161] about the possibility to realize a tt— junction for some 
values of parameters such as h, Q etc is not justified. 

Fig.l shows the contributions of the bulk singlet (/s) and Sz = triplet (/o) components to the critical current and 
corresponds to the tunnel Hamiltonian method. Fig. 2 displays the contributions of the singlet (/s), Sz = triplet 
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FIG. 3: The total critical current versus the exchange field obtained on the basis of the bulk Green's functions (solid line) 
and of the Green's functions taken at the MS/I interface (dotted line). The solid line corresponds to the tunnel Hamiltonian 
method. The normalization units are the same as in Fig.l. 



(/o) and \Sz\ = 1 triplet (Fi) components to the critical current. The solid and dotted curves in Fig. 2 are normalized 
partial critical currents defined as 

^3 = ^E/3(OMo = -^E/o(0)> ^i--^E^i(O) (29) 

UJ iU UJ 

where 7r^/8 — J^uj (2" + 1)^^ is the normahzation factor. The functions fi^oi^) ^^^^ ^i(O) are given by Eqs. (|25ll26)) . 
The lower (upper) dotted lines are due to the 3^ = and \Sz\ = 1 triplet components. It is seen that the current 
due to the singlet component decreases with increasing Xh, and the currents io,i due to the triplet components increase 
with increasing Xh- Interestingly, the current ii caused by the triplet component with nonzero projection of the total 
spin on the local z-axis is much larger than the current io caused by the Sz — triplet component. Meanwhile the 
current ii is absent in the tunnel Hamiltonian method at all (compare Figs.l and 2). 

In Fig. 3 we show the dependence of the total normalized critical current ic — is + "io + ii on A^j for = on the 
normalized exchange field Xh (dotted line). We compare this dependence with the dependence ic = ?3 + «o (solid line), 
i.e. with the critical current given by the tunnel Hamiltonian method, where 13.0 are determined by Ea. (|29p with 
/1 0(0) replaced by /|_q. One can see a significant difference between these dependencies. 



IV. CONCLUSIONS 



We have studied the influence of boundary effects on properties of magnetic superconductors with a spiral mag- 
netic structure. We used the well developed method of quasiclassical Green's functions. These functions obey the 
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Eilenberger (or Usadel) equations generalized to the case of an exchange field h acting on spins of free electrons and 
varying in space. For simplicity, we considered the case of temperatures close to the critical one, T^. Then, one can 
linearize equations for the condensate matrix Green's functions /. Due to a spatial dependence of the exchange field 
h, coefficients in the Eilenberger (Usadel) equations depend on the coordinate x. We excluded this dependence via a 
transformation which is equivalent to introducing a rotating coordinate system. In this local coordinate system the 
field h has only the z-component and does not depend on x. Solving these equations with corresponding boundary 
conditions, we have shown that near the boundary of MS with vacuum or an insulator, the condensate functions / 
differ essentially from their bulk values. 

In the rotating coordinate system, there are two components of the matrix /, /s and /o, in the bulk. These 
correspond to the singlet component and the triplet component with zero projection of the total spin on the z-axis. 
Due to boundary conditions, the corrections (5/o.3 to the bulk functions, /o,3, arise near the boundary, which are not 
small in comparison with /o,3. Besides, the triplet component Fi with nonzero projection of the total spin of Cooper 
pairs appears in the vicinity of the surface on the scale of the coherence length. The corrections (5/o,3 and function Fi 
oscillate with the period 271 /Q in space and decay inside the bulk over a length of the order of — v/^irT (ballistic 
case) or = \/ D/2ttT (diffusive case). The amplitude of the singlet component fj, decreases at the surface resulting 
in a suppression of the order parameter A near the surface. 

As an example of importance of the surface effects in MS's, we considered the dc Josephson effect in a MS/I/MS 
tunnel junction. The critical Josephson current Ic can be expressed in terms of components /o.3(0) and i^i(O) at the 
MS/I interface. The results are compared with the ones which are obtained on the basis of the tunnel Hamiltonian 
method and expressed in terms of the bulk condensate functions /o,3. This method was used in Ref.[16]. Although the 
formulae for Ic in Ref. ^IG'I and in this paper are similar, there is an essential difference between them. In the tunnel 
Hamiltonian method, the coefficient in front of cos 9 is the squared amplitude of the triplet Sz = Q component, /q . 
In fact, this coefficient is equal to /q (0) + -^i (0) (see Ea. (|28)) ). where /o(0) is the amplitude of the triplet component 
with zero projection of the spin on the local z-axis and Fi{Q) is the amplitude of the 15*2 1 — 1 triplet component at the 
interface. It turns out that, at least near Tc, the amplitude i^i(O) is much larger than /o(0). The tunnel Hamiltonian 
method can be applied to MSs only if the wave vector of the spiral, Q, is small enough: vQ « h (ballistic case) 
or DQ^ « h (diffusive case). However in this case the exchange field h should be small: < A (T << A) or 
h « [Tc — T)/T (A << T). Otherwise superconductivity will be destroyed. In this limit of small Q, the junction 
MS /I /MS is eq uivalent to the FS/I/FS junction. The Josephson current in FS/I/FS junctions was calculated in 
Refs.|33,[3l|,[3l|33. 

The surface effects may also change other characteristics of MS's such as the density-of-statcs (DOS) etc. Our 
consideration is restricted with temperatures T near Tc, where the DOS is close to that in the normal state and the 
variation of the DOS due to the surface effects is small. The calculation of the Green's functions in a finite system at 
low T is a more complicated task because the corresponding equations, strictly speaking, can not be linearized. This 
problem is beyond the scope of this paper. 
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VI. APPENDIX 

Here we obtain a formula for the Josephson current Ij in a MS/I/MS tunnel junction. We consider magnetic 
superconductors MS with a spiral magnetization described by the angle a{x) — Qx + 9 (right superconductor) and 
a{x) = Qx (left superconductor) so that 6 is the angle between the magnetization vectors at the MS/I interface. In 
order to obtain the expression for /j, we employ the boundary conditions [20. ,34. ,35.] 

fldjl = {2<jRB)-^[flJr] (30) 

where fi^r are the condensate functions in the left (right) superconductor, a is the conductivity of the superconductors 
in the normal state, and Rb is the junction resistance per unit area. The superconductors are assumed to be identical. 
The current is equal to [l^] 



10 



/ - (5a/8)^(27rT) ^Tr{f3 (g> ao (g> fi (g> dji} = ^^i(2^T) ^ Tr{f3 ® ao ® [/,,/.]} (31) 

UJ UJ 

where to = 7rT(2n + 1) is the Matsubara frequency and all the functions are taken at the interface {x = 0). 

We assume that the phase of the left superconductor is (p and the phase of the right superconductor is zero. Then, 
we can express the functions fi^r iu terms of the functions / found above with the help of transformations 



// =^ [/^0J7,0//®C/+®f/+, fr =^ Ur®fr®U+ (32) 

Here tJ^p = cos(</j/2) + i% (g) ao sin(y)/2) is the transformation matrix which relates a state with phase equal to zero 
and a state with a finite phase ip [l^]; Ui,r = cos(a;^r/2) + it^ ® cri sin(a;^r/2) with ai = Qx and = Qx + 6. Then, 
we substitute expressions ([5^ together with p?)) into Eg. ipTI) . Calculating the commutator in Eq. (PT|) . we come to 
Eq.lUS]). 

It is worth noting that the tunnel Hamiltonian leads to the same formula as Eg. pip if the functions fr^i are replaced 
by the bulk solutions, /o,3. 
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